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Abstract 



^^ We present results on the non-linear dynamics of inhomogeneous cosmo- 

logical models with irrotational dust and a positive cosmological constant, 
considering, in particular, a wide class with vanishing magnetic Weyl tensor. 
For those patches of the universe that do not recollapse to a singularity we 
QQ find a unique attractor, representing a de Sitter vacuum phase. For the (re- 

^~~^ )collapsing regions we find a family of (Kasner) attractors, so that generically 

these regions fall in spindle-like singularities. These results give substantial 
support to the idea that the universe can be very inhomogeneous on ultra- 
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in 



(^ large, super-horizon scales, whit observers living in those (almost) isotropic 

\.^ regions that emerge from an infiationary phase. 
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Cosmology in the XX century has been essentially based on the 2 (+1) parameter (i/o, 
Oo, and go) standard Friedman-Robertson- Walker (FRW) models. In 1981 inflation came 
on the scene as a possible solution to the conundrums of the big-bang FRW scenario [1,2]. 
However, despite the fact that these flaws are in a way or another related to the question 
FRW models cannot answer: why the observable universe looks isotropic^ and the inflationary 
scenario was proposed also to answer this question, in practice most of the work on inflation 
has been done in the framework of isotropic FRW models. It was indeed soon recognised 
that during an inflationary phase initially present small perturbations are swept away, a fact 
that led to the conjecture that this could be the signature for a more general property of 
inflation, going under the name of cosmic no-hair theorem [3]. Roughly speaking, inflation 
should erase previously present inhomogeneities, leaving us with a unique possible observable 
universe: the isotropic one. Investigations to prove some restricted version of this general 
conjecture were done first in the framework of homogeneous anisotropic models, [4] then also 
considering inhomogeneous spacetimes, [5] but with practical examples limited to geometries 
with some degree of symmetry (e.g. [6,2] and references therein). The dominant perspective 
emerging from these analysis is that in a universe model either infiation occurs everywhere, 
or there is no infiation at all. Instead, with truly inhomogeneous initial conditions, one can 
expect that there will always be patches of the universe that will infiate and isotropize, while 
others will not, eventually recoUapsing to a singularity. In this perspective, we propose that 
a weak cosmic no-hair theorem holds. 

It is the aim of this Letter [7] to present some results about the evolution of inhomoge- 
neous universes with irrotational dust of density p and a positive cosmological constant A, 
in particular, studying the dynamics of the wide class of models with vanishing magnetic 
part of the Weyl tensor, Hab = 0. 

These spacetimes were first considered in [8], while their first cosmological implementa- 
tion was given in [9], where it was also shown that for the case of dust their time evolution 
is given by a system of six first-order ordinary differential equation for p (matter density), 
(expansion scalar), CTi and (T2 (two independent eigenvalues of the traceless shear tensor crab)^ 
and El and E2 (eigenvalues of the traceless electric Weyl tensor Eab)- Thus, the evolution 
of each "fluid element" in the Hab = models is purely local once initial condition satisfying 
the appropriate constraint equations [8] are given. Because of this property we dubbed these 
models silent universes [10,11]. As shown in [10], Hab is a good approximation (at least at 
second order in perturbations of a FRW background) outside the Hubble horizon, where 
also pressure gradients can be neglected, so that these models may provide a fair picture of 
the universe on ultra-large scales. In other words, the time evolution of each super-horizon 
sized "volume element" of the universe is well described by the equations of the Hab = 
models. 

The dynamics of these models with A = has been presented in detail elsewhere [11]. 
Here we only point out that they put the flatness problem [1,2,12] in a rather unusual 
perspective: indeed, independently of its initial value, ultimately tends to zero, both for 
expansion and collapse, with most of the universe volume dominated by expanding voids. 
Deflning a± = |(cri it (T2) and E± = h{Ei it E2) the dynamics of irrotational dust with 
Hab = and A > is given by 
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where cr^ = 3(7+^ + (7_^ is the shear magnitude. Setting (7_ = £"_ = one obtains the 
dynamics of a set of universe models generalizing those of Szekeres to the case of A > 0; 
some of these models were explicitly given in [13]. [14] Equations (la) and (lb) are actually 
more general than for the case Hah = 0, as they hold for irrotational dust in general. Also, 
in this case the curvature of the 3-surfaces orthogonal to the matter 4-velocity [15] u" is 



(3)i? 



--02 + 2a2 + 2p + 2A. 
o 



Proofs of cosmic no-hair theorems for homogeneous [4] and inhomogeneous [5] spacetimes 
were essentially based on the assumption ^^'R < 0, although inflation has been found also 
in some ^^^ R > models [6,2,16]. For ^^^R < 0, it follows that > VSA, while < 
A — |0^ < 0, so that if at t* (an arbitrary initial time) 0* > (the model is initially 
expanding), then the model expands forever, with squeezed between the lower and upper 
bound and — > v 3A exponentially, thus approaching (at least locally) a de Sitter spacetime 
in a timescale a = a/3/A. 

Let us now drop the assumption ^"^'R < and make some general (i.e. not restricted to 
the Hab = case) remarks. From (lb) and 0* > v3A the upper bound 
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still holds, as in [4,5], but for 0* < v3A this becomes 
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Defining as usual [17,11] by 3i/i = the local scale factor £, from (la) one has that 
/9 — > as long as the expansion proceeds, and in this case from (lb) and (3), (4) one gets 
that — > V 3A if also anisotropy is suppressed, i.e. cr — > 0. However, is no more bounded 
from below in the general case, so that in general one can also expect recoUapse, even after 
an inflationary phase. From (lb) and (2) we see that if A > |0^ then ^"^'R > 0, and also 
if > then ^^^R > 0, while if A < ^0^ then < 0. In particular, this imphes that 



< — v3A is a no-return region: every trajectory entering it will undergo collapse with 
^"^'R > sufficiently close to the singularity. But if a patch of the universe starting with 
< V 3A has to inflate and isotropize, then in order to get the — > v 3A asymptotic value 
it has to super-inflate [f8], with > and ^^' R > 0. 

In terms of the usual density parameter = Om + f^y, with the matter and vacuum 
density parameters defined as usual, Om = 3/9/0^ and Oy = 3A/0^, it follows from (2) that 
> 1 ^ (3)i? > and ^^'^R < ^ < 1, and from (lb) that > ^ > 1, but the 
reverse relations do not hold, because of the shear. 

Now, let us consider the specific case of the dynamics of the Hab = models, described 
by system (1). For A = this system has a single stationary point, [19] given by the origin in 
phase-space, p = Q = cr± = e± = 0^ and corresponding to Minkowski spacetime. For A > 
this point bifurcates in 10 new points. Of these, 3 are static (0 = 0), with only one physically 
meaningful (the other 2 have /9 < 0), and represents the Finstein universe, with p = 2A and 
a± = E± = 0. Another point represents the spatially fiat de Sitter spacetime, with = v 3A 
(assuming expansion), and p = cr± = E± = 0. Two other points have (T_ = £"_ = 0, thus 
are degenerate and represent an oblate and a prolate configuration expanding exponentially 
with = a/A/3 and = v A respectively. Finally, the two other pairs of stationary points 
are physically equivalent replicas of these two latter. Now, the linearized stability analysis 
of these stationary points shows that all of them are saddle points, thus unstable, except 
the one representing de Sitter, which is asymptotically stable if > 0, i.e. de Sitter is the 
unique attractor for expansion. It then follows that those patches of the Hab = models 
that do not recoUapse expand toward a de Sitter phase and isotropize, i.e. cr-t — > and 
E± — > 0, with — > 1 and ^"^'R — > 0. Therefore, even in case the portion of the universe 
undergoing infiation was initially negligibly small, in a timescale of order a = a/3/A most 
of the universe volume will be isotropized by a phase of de Sitter-like infiation: a picture 
which is in the spirit of chaotic infiation [1]. According to the infiationary scenario, this de 
Sitter period will be ended by a suitable reheating process, after which a standard FRW 
phase will occur with close to unity. 

Another view of the phase-space for the Hab models is achieved using Om and Oy, 
together with the other dimensionless variables [20] Ej- = cr±/Q and e± = E±/Q'^. Obviously 
these variables diverge for = 0, so that one has to separately consider the > and 
< cases. Also, it is convenient to introduce a new "time" r = ib31n£, using the minus 
(plus) for the < (0 > 0) case, so that dr/dt > in both cases. 

Denoting by a prime the derivative with respect to r, by E^ = 3 E_|_^-|- E_^ the magnitude 
of the dimensionless shear, and by fto = Om — 20y the effective gravitational mass density 
parameter, the evolution equations for our variables for < read 

0' = |(2 + 12E2 + Og), (5) 
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OV = -^(2 + 12E2 + Og), (6a) 

f^M = -%^(12S2-l + Oa), (6b) 



SV = -i( 2 - 12 E^ - 6 E+ - Og) + -S_' + e+ , (6c) 

D O 

E'_ = ^(2- 12E2 + 12E+ -OG) + e- , (6d) 

e'+ = ^( 1 - 12 E^ + 9 E+ - Og) - S_ e_ + iE+ Om , (6e) 

6 D 

e'_ = ^( 1 - 12 E^ - 9 E+ - Og) - 3 E_ e+ + -E_ Om • (6f) 
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We see that now the Raychaudhuri equation (5) is decoupled from the rest of the system (6), 
i.e. all other equations do not depend on 0. The advantage we get in introducing system (6) 
is that, in addition to the non-static (0 ^ 0) stationary points of system (1), we also have 
a set of stationary points representing models for which the variables in (1) diverge , while 
those in (6) obviously have finite constant values. More precisely, system (6) admits many 
isolated stationary points, 3 of which are unphysical (Om = —3), with many degenerate 
physically equivalent triplets, and 2 physically equivalent sets of points parametrized by E_|_. 
The physically significant and distinct points are listed in Table I; points LI, L II and L III 
have ftv ^ and are the same non-static points of system (1) described above. Point D I- 
D VI are degenerate points with Oy = (i.e. for finite they represent A = models), and 
the set T III represents triaxial configurations parametrized by E_|_ (see the notes in Table 

I)- 

The linearized stability analysis shows again that for expansion > point L I, locally 

representing a fiat de Sitter universe, is again the unique attractor, while for collapse, < 0, 

the set Till, together with its conjugate Till (given by the sign exchange E_ — > — E_ and 

e_ — > — e_ in the expressions in Table I), is attracting. If collapse occurs — > — oo and 

Oy — > 0, while E-t and e± tend to finite values, and Om — > as well: it can be shown [11] 

that the sets Till and Till are locally equivalent to the Kasner models, thus the outcome 

of the stability analysis for system (6) is that collapsing configurations generically tend to 

a triaxial Kasner-like spindle singularity, with matter and A having no effects in the final 

stage. 



TABLES 

TABLE L Stationary points of system (6) (only the physically interesting points are listed). 
Eor expansion (0 > 0) only the de Sitter point L I is asymptotically stable, and for collapse (0 < 0) 
the only stable points are those of the Kasner set, i.e. the Till family and its conjugate THE All 
other points are saddles. Models marked with M are equivalent to Minkowski. 
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To get a clue of what goes on locally in an inhomogeneous universe with an effective 
cosmological constant we plot in Eig. 1 five significant cases of several numerical integrations 
of system (1). In Eig. 1 (left-up) we plot 0/v 3A, the upper bound (3) locally representing 
the open de Sitter universe, the = v 3A line representing the fiat de Sitter, and the bound 
(4) representing a closed de Sitter. In the strip |0| < v3A ^^'R > 0. The dynamically 
normalized, dimensionless 3-curvature scalar ^^'R/Q'^ is plotted in the right-upper corner. 
In the lower part we plot = Om + ^v and the deceleration parameter q = —30/0^ — 1: 
one has g < during infiation, q < —1 during super-infiation, and g — > 2 during collapse. 
When the de Sitter solution is asymptotically approached — > 1 and g — > — 1. We see 
that by and large for infiationary evolutions 0/v 3A gets below the fiat de Sitter line, then 
approaching it from below, asymptotically going as (4), with vanishing positive 3-curvature, 
^^'R — > 0"*". Related to this, we have observed in our numerical results that ^^'R does not 
change sign during the evolution, although we cannot formulate this as a rule; also, note 
that although ^^'R — > oo in the collapse, ^^'R/Q'^ — > 0: the approach to the singularity is 
Weyl dominated. Also, we observe that in general does not approach unity even during 
an infiationary phase, unless super-infiation occurs. 
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FIG. 1. Evolution of live inflationary cases, with one flnally re-collapsing, against \/A/3{t-\-to), 
where tg and the upper bounds (thick lines) can be read from (3) and (4). For the re-collapsing 
case, before to diverge at turn around (0 = 0), ^^'i?/©^ and Q are out of scale. 

In this Letter we analyzed Hab = <-^ab = p = models with a cosmological constant 
A > thought as representing the local value of the vacuum energy density, and possibly 
causing the occurrence of an inflationary phase in certain patches of the universe with 
suitable initial conditions. Of course in a realistic scenario such a cosmological constant 
will ultimately decay into radiation, reheating the universe to a standard ERW phase. The 
picture that emerges is that of a generally inhomogeneous universe, with large patches of 
it where, thanks to an early inflationary evolution, the local properties are close to that 
of a flat ERW model. In this picture, however, the observable parameters i^o, ^o and go 
should be thought as local values, not to be interpreted as giving the global properties of 
the universe. 
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